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ABSTRACT 


This  paper  presents  computer  estimations  of  the  response  of  models  of  mechanically  driven 
regularly  ribbed  fluid  loaded  panels.  The  estimates  are  presented  in  terms  of  the  spectral 
distribution  of  the  acceleration  of  the  panel  and  the  spectral  distribution  of  the  pressure  in  the  fluid. 
The  ribs  are  defined  in  terms  of  line  and  line  moment  impedances;  these  may  be  mass,  resistance, 
and/or  stiffness  controlled,  and  of  various  values  and  combinations.  Typical  estimations  of  the 
spectral  distribution  of  the  magnitudes  of  the  acceleration  are  displayed  as  functions  of  the 
normalized  wavenumber  that  lies  in  the  plane  of  the  panel  and  normal  to  the  ribs  and  the  normalized 
frequency.  The  magnitudes  of  the  acceleration  vanish  at  the  sonic  loci  when  fluid  loading  is 
included.  A  sonic  gorge,  with  a  nadir  at  the  sonic  loci,  characterizes  the  acceleration  of  the  fluid 
loaded  panel,  whether  unribbed  or  ribbed.  If  the  magnitudes  of  the  line  impedances  (and,  when 
appropriate,  also  the  line  moment  impedances)  of  the  ribs  are  not  unusually  high,  the  magnitudes 
of  the  acceleration  and  their  patterns  elsewhere  in  the  spectral  domain  remain  substantially  similar 
to  ‘hose  obtained  in  the  absence  of  fluid  loading. ^The  dissimilarities  are  readily  accounted  for  by 
straightforward  arguments  relating  to  the  increase  in  the  surface  mass  and  the  radiation  damping 
that  manifest  fluid  loading  on  panels.  The  corresponding  magnitudes  of  the  pressure  on  the 
surface  of  the  panel  do  not  vanish  at  the  sonic  loci  as  do  the  magnitudes  of  the  acceleration. 

Rather,  a  moderate  sonic  peak  appears  at  the  sonic  loci  and  a  sonic  ridge  replaces  the  sonic  gorge. 
The  loss  of  subsonic  components  in  the  pressure  on  a  plane,  as  the  plane  is  removed  from  the 
surface  of  the  panel,  is  clearly  displayed.  ^ 
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INTRODUCTION 


A  recent  paper  formulated  the  acoustic  behavior  of  ribbed  panels  immersed  in  uniform 
environments  and  presents  computer  estimations  that  examined  and  illustrated  some  aspects  of  the 
formalism  [1].  The  response  of  a  ribbed  panel  to  a  specific  external  drive  is  presented  in  terms  of 
the  spectral  distribution  of  the  velocity  (acceleration  or  displacement)  on  the  panel  and  the  pressure 
on  a  plane  at  or  above  the  panel.  An  external  drive  is  a  drive  which  is  independent  of  the  response 
of  the  dynamic  system  that  it  excites.  In  this  paper,  the  specific  external  drive  is  a  line  drive  that  it 
is  oriented  parallel  to  the  ribs.  The  estimations  that  are  computed  can  be  employed,  for  example,  in 
the  design  of  spectral  (wavevector-frequency)  filters  destined  to  measure,  in  real  situations,  similar 
responses  on  similar  dynamic  systems.  In  a  companion  paper  such  considerations  were  conducted 
on  panels  that  were  membranes  and  the  environmental  loadings  were  ignored  [2].  A  major 
emphasis  in  Reference  2  was  placed  on  the  comparison  between  panels  to  which  a  finite  and  an 
infinite  number  of  ribs  were  attached.  (It  would  help  the  reader  to  be  familiar  with  the  material 
discussed  in  References  1  and  2.)  This  paper  considers  a  more  elaborate  model  of  a  ribbed  panel 
which  may  be  plate-like,  and  fluid  loading  is  included.  However,  only  regularly  ribbed  panels  are 
considered;  i.e.,  panels  to  which  an  infinite  number  of  identical  and  regularly  spaced  ribs  are 
attached.  For  a  regularly  ribbed  panel,  one  derives  the  spectral  distribution  of  the  velcv  ty 
V  (k,  (D2)  terms  °f  tiie  impulse  response  function  G  (k  |  k ,  co2)  and  the  external  drive  Pe(k ,  o)2) 
both  expressed  in  the  spectral  domain;  namely, 


V(k,  ft>2)  =  J*G(k  |  k',  fi)2)  dk'  Pe(k',  »2);  A(k,  ©2)  =  icoV(k,  jo2); 
D(k,  ©2)  =  (ito) 1  V(k,  co2), 


G(k  |  k',  C)2)  =  GJk,  q2)  {5(k-k')  -  Ss( k  |  k',  @2)}  , 


(2) 
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Ss  (k  |  k,  ©2  )  =  [1  +  X  fL  (k  +  Kjt  0)2)]  1  H^k ,  &2) 


J^S(k  +  Ka-k') 

n 

(Kj)'1  X  exp  [  ixn  ( k  -  k')  ] 


.(3) 


[cf.  Eqs.  (8),  (9),  and  (lib)  of  Reference  2.]  In  Eq.  (1)  the  spectral  distributions  of  the  velocity 
V  (k,  g>  2),  the  acceleration  A  (k,  0)  2) ,  and  the  displacement  D  (k,  co  2)  on  the  panel  are  related. 
The  relationship  is  simple  so  that  using  one  form  of  the  response  of  the  panel  rather  than  another  is 
a  matter  of  convenience.  In  Eqs.  (2)  and  (3),  Goo(k,  ©2 )  is  the  surface  admittance  of  the  unribbed 
panel,  xn  is  the  position  in  the  x-domain  of  the  (n)th  rib,  K-  is  the  separation  wavenumber,  Kj  = 
jKj,  Kj  =  (2it/b),  where  b  is  the  separation  between  two  adjacent  ribs,  and  Hto  (k,  cg2)  is  the 
product  of  the  equivalent  surface  impedance  of  a  rib  and  the  surface  admittance  of  the  unribbed 
panel.  (When  the  ribs  are  absent,  b  serves  merely  as  a  linear  spatial  scale  factor.)  The  explicit 
expression  for  (k,  ©2)  is 

Ho*  (k,  tp2)  =  (Kj/271)  [  Z  (@2)  -  ik  Z^)  ]  Gro  (kco2  )  , 

where  Z  (o^)  and  Zm(o)2)  are  the  line  impedance  and  the  line  moment  impedance  of  a  typical  rib. 
(The  line  moment  impedance  Zm(®2)  a  viable  quantity  only  when  the  panel  can  support  a 
moment  response;  e.g.,  a  plate  responding  in  flexure.  Otherwise,  Zm(@2)  needs  to  be  set  equal  to 
zero  [2].)  The  spectral  vector  o>  2  =  {ky,to}  designates  the  Fourier  conjugate  of  y  and  t;  y  is  the 
spatial  variable  that  lies  along  the  ribs  in  the  plane  of  the  panel  and  t  is  the  temporal  variable;  see 
Fig.  1.  The  dependence  of  a  quantity  on  the  spectral  vector  variable  Q  2  is  suppressed  at  times  as  a 
matter  of  abbreviation;  e.g.,  G^fk,  0)2)  =  G^fk) .  Note  that  Ss  ( k  |  k ,  ©2) ,  stated  in  Eq.  (3),  is 
aliased  in  the  wavenumber  k  with  respect  to  the  harmonics  of  the  separation  wavenumber  Kj.  This 
aliasing  is  transmitted  also  to  the  drive  in  lieu  of  the  ribs;  the  aliasing  in  the  latter  quantity  was 
discussed  briefly  [2].  The  aliasing  in  k  with  respect  to  the  harmonics  of  Kj  may  be  stated  in  the 
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form  Ss  ( k  |  k  )  s  Ss(  k  +  Kj  |  k' ) .  [cf.  Eq.  (13)  of  Reference  2.]  The  external  drive  employed  in 
this  paper  is  of  the  form 


Pe  (k'>  ®2>  =  Pea  (®2>  eXP  (ik'xa  ) 


(5a) 


This  form  of  the  external  line  drive  qualifies  as  follows:  The  external  line  drive  pe(  x,  cp  2  )  is  assumed 
to  be  applied  centrally  at  xa  so  that  about  this  position  its  description  is  p«;(  x  -  xa,  £>2  )  •  It  follows 
that  pe(  x.  co  2  ) s  Pce(  X  -  xa,  £)  2  ).  Then,  by  a  Fourier  transformation,  one  obtains 
Pe  (k,  <o2)  =  Pee(k,  @2 )  exp  (ikxa),  where  typically  one  defines  Pee(k)  =  (2n)'u2  j  dx  pce(x)  exp  (ikx). 
Furthermore,  if  the  external  line  drive  is  localized  at  xa  so  that 


PeeC  X  -  Xa,  @  2  )  =  (27t)  ^  Pea(®  2>  5(  X  "  Xa)  ’ 


(5b) 


then  Pee(k,w  2)  =  Pea({02)  and  Eq.  (5a)  is  explained.  From  Eqs.  (1)  through  (3)  and  (5)  one 
obtains 


V(k,co2)  =  Voo(k,(g2){  1- 


[1+2  +  Kj,  @2  )]  1  £  Hoo0c  +  ^  ©2  )  expo^)}  » 


where 


VJk,  @2 )  =  co  2 )  Peatoj  )  exp(ixak) 


(6) 


(7) 


PARAMETRIC  DEFINITION  OF  THE  FLUID,  THE  PANEL,  AND  THE  RIBS 


To  facilitate  the  estimation  of  V(k,  co  2)  [A  (k,  co  2)  or  D  (k,  £>  2)  ]  as  stated  in  Eqs.  (1)  and 
(6),  it  remains  to  express  more  explicitly  the  quantities  and  parameters  that  are  involved  in  this 
equation.  The  surface  admittance  of  the  panel  is  expressed  in  the  form 
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G^ck,  (j^)  =  [Zp(k,  o2  )  +  Z/k,  @2)]  1  ;  Z/k,  ©2)  =  Z^(k,  co2)  +  Z^/k,  co2), 


(8) 


where  Z/k,  0)2)  is  the  loading  on  the  panel  by  the  uniform  environment ,  Z^/k,  g^)  is  this 
loading  on  one  (top)  surface  and  Z  * /(k,  g^)  is  this  loading  on  the  other  (bottom)  surface  of  the 
panel,  and  Zp(k,  ©2 )  is  the  mechanical  surface  impedance  of  the  panel  [1]. 

In  this  paper  the  environmental  loading  is  limited  to  that  of  a  semi-infinite  fluid  occupying 
the  space  above  the  panel  as  shown  in  Fig.  1.  The  fluid  loading  on  the  panel  may  then  be 
expressed  in  the  form  [1,2] 

Z  Jk,  g),)  =  Zt/k,  co,)  =  pco/k3  =  pc  (k3c/ co)*1  ,  (9) 


where  p  is  the  density  and  c  is  the  speed  of  sound  in  the  fluid  and 

-  1/2  ~  -  1/2  T 

(k3c/co)  =  [  1  -  k2  ]  u  [  1  -  r  ]  -  i  [  k2  -  1  ]  U  [  k2  -  1  ] 
ti2  =  ( c  /  co  )2  ( k2  +  k^ ) 


GO) 


In  this  paper  the  mechanical  surface  impedance  of  the  panel  is  stated  in  the  simple  but 
orthotropic  form 


Zp(k,®2)=  icom  [  1  -  { (k  /  kp  )2  +(ky  /  kpy  )2  }p  ]  , 


(11) 


where  m  is  the  mass  per  unit  area  and  {  kp,  kpy }  is  the  free  wavevector  on  the  panel.  The  panel  is 
membrane-like  if  p  =  1,  and  is  plate-like  responding  in  flexure  if  p  =  2.  The  free  orthotropic 
wavevector  kp  is  written  in  the  form 

kp  ={kp,kpy}  ;  kp  =  kpo  (1  -  iT|p  ) ;  kpy  =  kpyo(l  ■  i”Hpy  )  »•  Kpo  ={kpo,kpyo}  t  (12) 


where  {t)p,  Tipy}  is  the  loss  factor  "vector”.  For  a  plate  and  a  membrane  that  is  made  to  simulate  a 
plate  responding  in  flexure 


5 


kpo  —  (CDOJc  /  c2);  l^yo  —  (OJCOcy  /  c2  ) ;  |kpol  —  (kp0  +  kpyo) 


(13) 


where  coc  and  G)cy  are  the  critical  frequencies  in  the  x-  and  y-domain,  respectively.  These  critical 
frequencies  are  defined  with  respect  to  the  speed  of  sound  c  in  the  fluid  facing  the  panel.  For  a 
plate  in  longitudinal  response  and  a  membrane  in  flexural  response 


kpo  —  I  coc  )2  (kc)2  (c  /  c^)2;  kpyo  —  (to  /  cocy)2  (kcy)^  (c  /  c  ^y)2> 


(14) 


where  {kc,  key}  =  {(ti)c  /  c),  (C0cy  /  c)}  and  the  free  wave  velocity  in  the  panel  is  {c^,  c^y}.  The 
surface  admittance  may  then  be  derived  from  Eqs.  (8),  (9),  and  (1 1)  in  the  explicit  form 


Gw  (k,  to y)  —  [Z^  ( k,  Q)  2)  ]  I 


2^P 


Z^  (k,  co2)  =  (i  (om)  [l-iOc/kpf^dCy/kpyfr-aE^  O^c/co,)'1]  , 


(15) 


where  the  fluid  loading  parameter  ^  is  defined 

Ec  =  (  p  c  /  (0cm ) 


(16) 


In  this  paper  the  panel  is  assumed  to  be  isotropic;  kp0  =  kpy0;  Tip  =  ppy,  o\.  =  tocy.  The  orthotropic 
form  for  the  mechanical  surface  impedance  of  the  panel  is  stated  to  indicate  that  simple  orthotropy 
can  be  readily  accommodated  in  the  estimations  of  the  response. 

The  line  impedance  and  line  moment  impedance  of  a  typical  rib  may  be  cast  in  various 
forms  some  of  which  may  be  elaborated  to  reflect  the  various  properties  that  the  ribs  and  their 
attachment  to  the  panel  may  possess  [1].  In  this  paper  the  explicit  expressions  for  the  line 
impedance  and  line  moment  impedance  are  cast  in  the  simple  forms 


Zj  =  Z=  i  0)  (  Mg+  Mo ) ;  Mo  =  M{l-  ( co/tUo  )2  (  1  -  i  r\0 )  }*J  , 


(17) 
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z Mj  -  i  ®  (M  /  kMo)  ;  kMo  -  kM  { 1  -  ( co  /  coMo  )2  (  1  -  i  TlMo  )  }  ( M  /  M0  )  , 

(18) 

respectively,  where  Mg  and  M  are  masses  per  unit  length,  kM  is  a  specified  wavenumber,  a>0  and 
como  316  given  resonance  frequencies,  and  and  "Hmo  are  assigned  loss  factors.  The  line 
impedance  of  the  rib,  stated  in  Eq.  (17),  is  slightly  more  elaborate  than  that  assumed  in  Reference  2 
in  that  a  mass  controlled  impedance  (icoMg)  is  added  and  placed  adjacent  to  the  panel  [1].  Also,  the 
line  moment  impedance  of  a  typical  rib  is  included  here. 

FromEqs.  (4),  (15),  (17),  and  (18)  one  obtains 

H«,  ( k.  a  2  )  =  [(Me/mb)+  ( M„/mb){  1  -  i(  k  /  kMo )}] 

[  l  -  {  (  k  /  ^  (  ky  /  fcpy  ^  lP  -  i^^C/O)/']  .  (19) 

Again,  the  elaborations  introduced  in  Eqs.  (17)  and  (18)  and,  in  turn,  in  Eq.  (19),  are  not 
employed  directly  in  this  paper;  they  are  included  in  order  to  indicate  the  readiness  with  which  they 
can  be  accommodated  in  the  estimations. 


ESTIMATIONS  OF  THE  RESPONSE  OF  AND  THE  PRESSURE  ABOVE  THE  PANEL 

Uniform  environmental  loadings  and  fluid  loading  in  particular,  can  be  readily  accommodated  in 
the  computations;  i.e.,  there  is  no  difficulty  inserting  Eqs.  (4)  and  (8)  in  general,  and  Eqs.  (15)  and  (19) 
in  particular,  in  Eq.  (6).  The  resulting  equation  evaluates  the  response  of  the  panel  merely  in  terms  of 
summations  over  harmonics  with  respect  to  the  separation  wavenumber  iq;  Kj  =  (2rc/b).  For  reasonable 
panels,  ribs,  and  fluid  loading,  the  summations  can  be  approximated  by  a  small  number  of  terms, 
therefore  the  computations  can  be  handled  by  a  desk  top  computer,  [cf.  Appendix  A.] 
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The  pressure  spectral  distribution  P(  k,  ©  2,  z)  in  the  fluid  on  a  plane  a  distance  z  above  the  panel 
is  derived  directly  from  the  velocity  spectral  distribution  V(  k,  ©2)  on  the  surface  of  the  panel;  namely, 


P(  k,  ©2,  z  )  =  exp  (-i  k3  z  )  Zf<  k,  fi? 2)  V(  k,  ©2), 

where 

Zf(k,  co  ,)  =  Zt/k,  co ,)  =  (pco  /  k3)  =  pc  (k3  c/co)'1  , 

and  use  is  made  of  Eq.  (9)  [1].  The  explicit  expression  for  (k3  c/co)  is  given  in  Eq.  (10).  Thus,  once  the 
velocity  spectral  distribution  of  the  panel  is  estimated,  the  pressure  spectral  distribution  on  a  plane  a 
distance  z  above  the  panel  readily  follows,  as  «tated  in  Eq.  (20),  where  Zf(k,  ©2 )  represents  the 
conversion  factor  from  the  spectral  distribution  of  the  velocity  V(k,  ©2  )  °f  the  panel  to  the  spectral 
distribution  of  the  pressure  P(  k,  ©  2,  0  )  [s  P(  k,  ©  2)]  in  the  fluid  on  the  surface  of  the  panel.  With  this 
interpretation,  the  factor  exp  (-ik3  z )  represents  the  filtering  performed  by  propagation  through  the  slab 
of  fluid  between  the  surface  of  the  panel  and  a  control  plane  above  it  where  the  pressure  is  estimated  f  1]. 
Eq.  (10)  shows  that,  for  a  supersonic  component  in  P  (k,  ©2),  this  filtering  is  a  mere  phase  change  — 
propagation  without  decay.  For  a  subsonic  component  in  P  (k,  ©  2),  this  filtering  is  an  exponential 
decay  with  an  exponent  that  increases  linearly  with  the  normal  distance  of  the  control  plane  above  the 
panel. 


(20) 

(21) 


NORMALIZATIONS  AND  DISPLAYS 

The  preceding  formalism  is  used  to  calculate  the  spectral  distributions  of  the  response  of  several 
representative  cases  of  interest.  The  response  quantities  chosen  are  the  acceleration  A(k,  ox, )  on  the 
surface  of  the  plate  [which  is  simply  related  to  the  velocity  V(k,  ©2  )  and/or  displacement  D(k,  ©,  );  see 
Eq.  (1)]  and  the  pressure  P(  k,  ©  2,  z )  in  the  fluid  on  a  control  plane  placed  a  normal  distance  z  above 
the  surface  of  the  plate.  The  relevant  spectral  distributions  are  cast  in  the  normalized  forms: 
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A  (k,  ©2)  =  A(k,  ft)2)  [Pea(jo2)/m]  1  =  (co/coc)  V(k,  ©2)  ; 

v(k,  ©2)  =  v(k,  «2)  rpca(e2)/i03cm] 1  * 


(22) 


P(k,  co2,  z)  =  P(k,  (o2,  z)  [Pea  (co  2)  (ec/i)]  ,  (23) 

where  e ■<.  is  defined  in  Eq.  (16).  Only  magnitudes  of  these  quantities  are  displayed.  By  and  large,  these 
magnitudes  are  displayed  clipped  by  certain  predetermined  and  convenient  clipping  values  (thresholds) 
so  that  only  magnitudes  that  exceed  these  thresholds  are  shown.  The  clipping  values  are  designated  by  a 
superscript  of  zero;  e.g.,  the  clipping  values  of  |  Aj  and  |P|  are  thus  designated  |A°|  and  |P° |,  respectively 
[2].  The  estimations  may  be  conveniently  displayed  on  planes  defined  by  any  two  of  the  normalized 
dependent  variables  (k/iq),  (ky/Kj),  and  (co/coc),  with  the  third  held  fixed.  In  this  paper,  the  {(k/Kj), 
(co/coc)}-plane  is  used  and  the  magnitudes  of  a  quantity  are  computed  as  a  function  of  the  normalized 
wavenumber  (k/Kj),  at  specific,  equal,  and  successive  values  of  the  normalized  frequency  (oVcoc)  with  a 
fixed  value  of  (ky/Kj)  [2].  Also,  the  computations  are  focused  on  limited  portions  of  the  {(k/Kj), 
(a>/(Dc)}-plane;  usually  the  lower  wavenumber  and  frequency  portions.  [The  graphs  are  displayed  in  the 
format  of  Figs.  3  and  4  of  Reference  2;  acquaintance  with  the  format  of  these  figures  will  assist  readers 
unacquainted  with  these  display  procedures.] 

COMPUTER  ESTD  ATIONS 

In  Section  13,  the  parameters  of  influence  on  the  estimated  magnitudes  of  interest  are: 
p,  (cA/c),  (ky/Kj),  (M/mb),  rip,(bkc),  (co0/  toc),  r\Q,  (bkM),  (coMo/  coc),  nMo,  (wc/tocy),  (c^  / c/y), 

(rip  /  rjpy),  (zb),  (xa/b),  ec,  |A°|,  [P  °|,  etc.  In  describing  the  figures,  it  is  convenient  to  define  a  set 
of  standard  conditions  for  the  estimations  of  the  spectral  distributions  of  the  magnitudes  of  the 
response.  The  standard  conditions  are: 
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The  panel  is  isotropic, 


(M/mb)  =  0.3,  (bkm)  =  106,  (bkc)  =  16,  (bky)  =  0  ,  (xa/b)  =  0.3  ; 

|A°!=  1.4,  |P°|  =  1.4,  and 

^  =  |5xl0-3  ,  p  =  1 

p  \2.5xl0-3  ,  p  =  2 


0.1  with  fluid  loading 

0.0  with  fluid  loading  removed 

(24) 


Note  that  the  standard  conditions  are  defined  with  the  thresholds  |A  °|and  |P  °|  equal.  The  standard 
conditions  are  conveniently  defined  so  that  only  departures  need  reporting  in  the  text  and  figure 
captions.  Figure  2  displays  a  set  of  computed  estimates  relating  to  the  acceleration  spectral 
distribution  of  line  driven  unribbed  panels.  In  Fig.  2a  the  acceleration  spectral  distribution  on  an 
unribbed  fluid  loaded  membrane  is  displayed;  the  membrane  is  made  to  simulate  the  dispersion 
properties  of  a  plate  responding  in  flexure,  [cf.  Fig.  4  of  Reference  2.]  The  expected  vanishing  of 
the  response  of  the  panel  at  the  sonic  loci  and  the  associated  sonic  gorge  is  clearly  discemable.  In 
Fig.  2b  the  corresponding  estimations  for  an  unribbed  fluid  loaded  plate  responding  in  flexure  are 
displayed.  The  analogies  and  differences  between  the  response  on  a  membrane  that  simulates  a  plate 
in  flexure  and  a  plate  in  flexure  can  be  deciphered  by  comparing  Fig.  2a  with  Fig.  2b.  Except  for 
minor  details,  the  patterns  in  the  response  show  that  a  membrane  can  be  analytically  modeled  to 
simulate  the  response  of  a  plate  in  flexure,  [cf.  Eqs.  (13)  and  (14).]  The  clipped  version  of  Fig.  2b 
is  displayed  in  Fig.  2c. 


The  corresponding  pressure  spectral  distributions  on  a  plate  responding  in  flexure  and  on  a 
plane  some  distance  above  it  are  shown  in  Figs.  3a  and  b,  respectively.  The  manner  in  which  the 
acceleration  spectral  distribution  converts  into  the  pressure  spectral  distribution  on  the  surface  of 
the  panel  is  illustrated  by  comparing  Figs.  2b  and  3a;  the  absence  of  the  sonic  gorge  in  the  spectral 
distribution  of  the  pressure  is  noted.  Indeed,  not  only  is  the  sonic  gorge  ironed  out  but  moderate 
peaks  occur  at  the  sonic  loci  and  a  sonic  ridge  is  formed.  The  suppression  of  subsonic 


10 


components  in  the  pressure  spectral  distribution  on  a  plane  parallel  to  the  panel  as  the  plane  is 
further  removed  from  the  panel  is  illustrated  by  comparing  Figs.  3a  and  b  [cf.  Eqs.  (10)  and  (20)]. 

The  computed  estimations  relating  to  the  response  of  a  regularly  ribbed  panel  that  is  mostly 
fluid  loaded  and  mostly  plate-like  is  now  introduced.  The  displays  are  in  the  format  of  Fig.  2c  in 
that  the  magnitudes  displayed  are  clipped.  In  Appendix  A,  the  formalism  stated  in  Eq.  (6)  admits 
readily  to  accounting  for  fluid  loading  and/or  to  accounting  for  regularly  ribbed  plates  responding 
in  flexure.  It  is  intended  just  to  illustrate  the  influence  of  fluid  loading  on  panels  that  are  plates 
responding  in  flexure;  other  kinds  of  panels  (e.g.,  membranes)  and  other  types  of  reponses  (e.g., 
longitudinal)  which  were  widely  used  in  Reference  2,  are  omitted  for  brevity. 

Figure  4  displays  the  clipped  magnitudes  of  the  acceleration  of  a  regularly  ribbed  fluid 
loaded  plate  responding  in  flexure.  In  this  Figure  the  line  impedance  of  each  of  the  identical  ribs  is 
assumed  to  be  mass  controlled,  and  the  line  moment  impedance  is  assun  to  be  negligible.  The 
clipping  value  in  all  of  Fig.  4  is  identical.  The  acceleration  spectral  distribution  of  a  standard  fluid 
loaded,  regularly  ribbed  and  isotropic  plate;  i.e.,  (M/mb)  =  0.3  ,  T]p  =  2.5  x  10~3,  (bkj  =  16,  coc  = 
cocy,  Tip  =  T]py,  and  £c=  0.1,  is  depicted  in  Fig.  4a.  Figure  4b  is  the  same  as  Fig.  4a  except  that 
fluid  loading  is  removed  by  setting  £c=  0.  [In  this  connection,  it  is  noted  that  a  fluid  loading 
represented  by  £c=  0.1  is  "heavy".]  The  similarities  between  the  displays  shown  in  Fig.  4a  and  b 
are  clear  enough.  However,  close  observation  shows  minor  dissimilarities;  e.g.,  the  dispersion 
peaks  occur  at  slightly  different  loci  in  Fig.  4a  and  b.  A  fairly  common  procedure  replaces  fluid 
loading  on  a  panel  by  adding  surface  mass  and  increasing  the  loss  factor  to  compensate  for  the 
radiation  damping.  How  adequate  is  such  a  procedure?  The  computation  performed  may  be  used 
to  answer  this  question  in  part.  A  surface  mass  of  (p/kpj)  is  added  to  the  surface  mass  of  the  panel 
and  the  distributed  damping  designated  by  T|p  is  increased  by  the  factor  y,  where  kpj  is  the  free 
wavenumber  on  the  panel  corrected  for  fluid  loading  and  y  >  1 .  To  accommodate  this  procedure, 
fluid  loading  is  ignored  and  the  mechanical  surface  impedance  of  the  isotropic  panel  is  replaced 
instead  by 
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Zp(k,  co2)  =  icom  [l  +  (p/mkpl)  -  {(k2  +  k^)/k2}P]  ; 
kp  =  kpo(1_i7T1P)  :  kJ  =  (cococ/c2)  - 


(25) 


where 


1  +  (p/mkpj)  =  (kp!^)215 


(26) 


[cf.  Eqs.  (1 1)  and  (15).]  If  one  finds  that  can  be  approximated  by  k^  in  Eq.  (26)  then  in  Eq. 
(25)  one  may  set 


(p/mkpj)  =  (p/mkp0)  =  (coc/co)1/2  ec  . 


(27) 


Figure.  4c  illustrates  the  use  of  this  procedure  for  the  standard  (moderate)  conditions  under  which 
Fig.  4a  is  obtained.  The  similarities  between  Figs.  4a  and  c  are  increased  and  the  dissimilarities  are 
decreased  as  compared  to  those  between  Figs.  4a  and  b;  e.g.,  the  dispersion  peaks  are  substantially 
coincident  in  Figs.  4a  and  c.  However,  closer  examination,  especially  in  the  unclipped  versions  of 
these  figures,  shows  that  the  sonic  gorge  is  not  present  in  Fig.  4c.  At  and  in  the  vicinity  of  the 
sonic  region,  the  similarities  between  Figs.  4b  and  c  are  closer  than  those  between  Figs.  4a  and  c. 
This  indicates  that  the  above  procedure  for  accounting  for  fluid  loading  may,  with  caution,  be 
useful  for  certain  purposes. 

In  Fig.  5  the  mass  ratio  (M/mb)  is  increased  from  0.3  to  the  rather  large  value  of  3: 
otherwise  the  standard  conditions  imposed  on  Fig.  4a  remain  intact.  In  Fig.  5a  the  fluid  loading 
remains  standard  at  £c=  0.1.  In  Fig.  5b  the  fluid  loading  is  removed;  £c=  0.  Comparing  Figs.  5a 
and  b  reveals  that  fluid  loading  does  significantly  influence  the  patterns  when  the  line  impedance  of 
a  typical  rib  is  "heavy."  This  kind  of  influence  of  fluid  loading  supports  the  notion  that  fluid 
loading  tends  to  "mend  and  soften"  the  discontinuities  associated  with  attaching  heavy  ribs  to  the 
panel  [3-5], 
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In  Fig.  6  the  line  impedance  of  a  typical  rib  is  decreased  and  the  line  moment  impedance  of 
a  typical  rib  is  introduced.  This  is  accomplished  by  choosing  a  mass  ratio  of  (M/mb)  =  0.05  and 
(kMb)  =  9.  [In  Figs.  4  and  5  the  line  moment  impedance  of  a  typical  rib  is  removed  by  setting 
(kMb)  =  106;  see  Eq.  (24).  ]  In  Fig.  6  the  line  moment  impedance  is  mass  controlled  as  is  the  line 
impedance.  In  Fig.  5a  fluid  loading  is  present  with  Ec=  0.1,  and  in  Fig.  5b  fluid  loading  is 
removed  by  setting  £c=  0.  The  similarities  and  dissimilarities  between  Figs.  6a  and  b  are 
analogous  to  those  found  between  Figs.  4a  and  b.  The  introduction  of  line  moment  impedances 
shows  no  dramatic  effects  in  the  response  of  the  panel  to  a  line  drive. 

Attention  is  now  focused  on  computer  estimations  of  the  pressure  spectral  distribution  on 
and  above  the  plane  of  the  regularly  ribbed  fluid  loaded  panel,  [cf.  Figs.  3a  and  b.]  Figs.  7a,  b, 
and  c  display  the  magnitudes  of  the  pressure  spectral  distribution  above  a  plate  responding  in 
flexure  and  under  standard  conditions;  see  Eq.  (24).  In  Fig.  7a  the  pressure  spectral  distribution  is 
that  on  the  surface  of  the  plate;  (zk<.)  =  0.  In  Figs.  7b  and  c  the  magnitudes  of  the  pressure  spectral 
distribution  are  those  on  control  planes  removed  by  (zlq.)  =  2  and  by  (zkj  =  12,  respectively, 
normally  off  the  plate.  The  suppression  of  the  subsonic  components  by  this  removal  is  clearly 
evident,  as  is  the  sonic  ridge.  Comparing  Figs.  3a  and  b  and  Fig.  7  illustrates  the  contribution  that 
ribs  make  to  the  pressure  spectral  distribution.  The  contribution  made  by  the  ribs  in  the  supersonic 
region  is  of  particular  interest  In  this  connection,  comparing  Figs.  4a  and  7a  again  reveals  the 
depression  at  and  in  the  vicinity  of  the  sonic  region  in  Fig.  4a,  and  the  enhancement  at  and  in  the 
vicinity  of  the  corresponding  spectral  region  in  Fig.  7a.  Indeed,  a  sonic  gorge  in  Fig.  4a  is  the 
counterpart  to  a  sonic  ridge  in  Fig.  7;  e.g.,  Fig.  7b  [6].1  The  dominant  influence  of  fluid  loading 

1  One  needs  to  differentiate  between  the  spectral  distribution  of  the  pressure  in  a  plane  above  the 
surface  of  a  vibrating  surface  and  the  pressure  radiated  to  a  localized  spatial  region  above  the 
surface.  These  two  manifestations  of  the  pressure  field  are  related  but  not  coincident  Indeed, 
characteristics  exhibited  in  one  may  not  be  exhibited  in  the  other.  Stretching  a  point  to  illustrate  the 
statement,  consider  a  panel  vibrating  under  a  local  external  drive.  The  sonic  ridge  will  be  present 
in  the  spectral  distribution  of  the  pressure  on  a  plane  removed  somewhat  from  the  panel,  [cf. 
Figure  2e.]  On  the  other  hand,  the  pressure  radiated  to  a  localized  region  on  that  plane,  but  that  is 
far  removed  from  the  location  of  the  external  drive  will  not  track  the  sonic  ridge,  notwithstanding 
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at  and  near  the  vicinity  of  the  sonic  range  is  thus  manifested;  the  influence  is  dominant  because  the 
fluid  surface  impedance  in  that  spectral  range  is  the  dominant  term  in  the  surface  impedance  of  the 
panel.  Again,  these  displays,  and  those  presented  in  Reference  2  were  given  to  briefly  establish 
the  kind,  the  type,  and  the  wide  range  of  computer  estimations  that  can  be  conducted  in  terms  of 
the  formalism  presented.  No  attempt  was  made  to  cover  a  specific  feature  of  interest.  This 
coverage  can  be  provided  when  a  specific  question  or  need  arises. 


that  the  contribution  to  this  radiated  pressure,  at  grazing  angles,  is  contributed  by  the  near  sonic 
components  on  the  plane  of  the  panel.  The  lack  of  tracking  is  explained  when  one  realizes  that  the 
radiated  pressure  is  contributed  by  the  appropriate  spectral  components  on  the  plane  of  the  panel 
that  are  related  to  the  velocity  and  not  to  the  pressure.  The  far  field  radiated  pressure  tracks,  by  and 
large,  the  spectral  distribution  of  the  velocity  and  not  of  the  pressure  [6]. 
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APPENDIX  A 


COMMENTS  ON  IRREGULAR  RIBBING 

It  is  crucial  to  realize  that  the  simplicity  in  the  computations  afforded  in  Eq.  (6)  is  in  direct 
consequence  of  the  regularity  of  the  ribs  and  their  attachment  to  the  panel.  For  irregularly  ribbed 
panels  this  simplicity  is  not  afforded.  If  the  panel  is  not  regularly  ribbed,  it  becomes  necessary  to 
evaluate  the  (transfer)  line  admittance  [1] 


goo  (x’  ©2)  =  J  Goo  (k.  © 2)  dk  exp  (-ikx) 


(Ala) 


Moreover,  when  the  panel  can  support  a  moment  response,  it  becomes  necessary  to  evaluate  also 
the  derivative  of  g^  (x,  50  2)  with  respect  to  x;  namely  [1] 

g'  (x,  to  2)  =  (3/3x)  g  (x,  co  2)  . 

00  (Alb) 


In  the  absence  of  fluid  loading  the  integration  involved  in  Eq.  (Al)  is  simple  enough,  but  is 
difficult  when  fluid  loading  is  included  [3-5].  Moreover,  since  the  response  of  the  ribbed  panel  is 
a  functional  of  g^  (x,  co2),  and,  when  appropriate  also  of  g^  (*,  ©2),  and  these  functions  must  be 
evaluated  over  an  extensive  number  of  values  of  x  which  are  not  simply  related,  the  evaluation  of 
the  response  is  beset  by  difficulties  that  are  compounded.  This  is  the  sense  in  which  Eq.  (6)  is 
considered  relatively  simple.  The  regularity  replaces  the  integral,  and,  when  appropriate,  its 
derivative  with  respect  to  x,  by  evaluations  at  a  number  of  discrete  values  of  the  wavenumber. 
Indeed,  one  may  claim  that  the  regularity  in  the  rib  spacing  not  only  facilitates  accounting  for  the 
fluid  loading  but  also  for  determining  the  response  of  a  plate  in  flexure.  In  the  case  of  a  plate 
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responding  in  flexure,  the  integral  in  Eq.  (Al)  is  more  complicated  than  it  is  for  a  mere  membrane. 
Finally,  the  regularity  of  the  separations  among  the  ribs  is  more  significant  than  the  identity  of  the 
line  impedances  and,  when  appropriate,  also  the  identity  of  the  line  moment  impedances  of  the 
ribs.  The  discrete  values  in  the  wavenumber,  which  are  introduced  by  the  harmonics  with  respect 
to  Kj,  are  not  destroyed  by  relaxing  the  identity  of  the  line  and  line  moment  impedances  of  the  ribs 
[2].  Thus,  it  is  conceivable  that  one  could  relax  the  requirement  for  the  identity  of  the  impedances 
of  the  ribs  without  entirely  forfeiting  the  advantages  that  an  equation  like  Eq.  (6)  holds  in 
accounting  for  the  fluid  loading  and,  when  appropriate,  in  accounting  also  for  "the  plate  in  flexure" 
properties  of  the  panel. 
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Fig.  1 .  A  sketch  of  a  regularly  ribbed  fluid  loaded  panel  showing  the  coordinate 
system  and  the  orientation  and  locations  of  the  ribs  and  the  external  line  drive. 


(k/Kj) 


Fig.  2.  Magnitudes  of  the  acceleration  spectral  distribution  for  an  unribbed 
fluid  loaded  panel  under  standard  conditions  [Eq.  (24)]  as  a  function  of  the 
normalized  wavenumber  (k/K})  at  successive  and  equal  increments  of  the 
normalized  frequency  (o)/(Oc). 
a.  A  membrane  that  simulates  a  plate;  p  =  1 . 
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(co/co  J)  =  0.6 


(0)/C0c)  =  0 


(k/K,) 


((0/(0c)  =  0.6 


((o/(oc)  =  0 


(k/K,) 


Fig.  3.  Magnitudes  of  the  pressure  spectral  distribution  on  a  plane  above 
an  unribbed  fluid  loaded  plate  under  standard  conditions  [Eq.  (24)]  as  a 
function  of  the  normalized  wavenumber  (k/K,)  at  successive  and  equal 
increments  of  the  normalized  frequency  ((o/coc). 

(a)  (k,)  =  0. 

(b)  (zkc)  =  12. 


Fig.  4.  Magnitudes  of  the  acceleration  spectral  distribution  of  a  regularly 
ribbed  plate  responding  in  flexure  under  standard  conditions  [Eq.  (24)]  as 
a  function  of  the  normalized  wavenumber  (k/Kj)  at  successive  and  equal 
increments  of  the  normalized  frequency  ((o/co^. 

(a)  £c  =  0.1. 

(b)  ec  =  0.0. 


Fig.  4  (continued) 

(c)  Introduction  of  fluid  loading  in  the  manner  specified  in  Eqs.  (25) 
through  (27)  with  e  =  0.1  and  y  =  2. 


V  — C' 


0 


b) 


3 


Fig.  5.  As  in  Fig.  4  except  that  the  mass  ratio  (M/mb)  is  increased  from  the 
standard  value  of  0.3  to  3. 


(a)  ec  =  0.1. 

(b)  ec  =  0.0. 
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Fig.  6.  As  in  Fig.  4  except  that  the  mass  ratio  (M/mb)  and  the  normalized 
moment  wavenumber  (bkM)  are  decreased  from  their  standard  values  of 

0.3  and  106,  respectively,  to  the  values  of  0.05  and  9,  respectively. 

(a)  ec  =  0. 1 . 

(b)  £c  =  0.0. 
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(co/coc)  =  0.6 


(co/co c)  =  0  - 


(o)/(o_)  =  0.6 


((0/co^  =  0 


(k/K,) 


Fig.  7.  Magnitudes  of  the  pressure  spectral  distribution  on  planes  above  the 
surface  of  a  regularly  ribbed  fluid  loaded  plate  responding  in  flexure  under 
the  standard  conditions  [Eq.  (24)],  as  functions  of  the  normalized  wavenumber 
(k/Kj)  at  successive  and  equal  increments  of  the  normalized  frequency  (co/coc). 
[cf.Fig.  3.] 

(a)  zkc  =  0. 

(b)  zk  =  2.  24 
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